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ABSTRACT 


Within  the  framework  of  the  "exact"  linear  theor-  an  important 
class  of  wave  propagation  problems  in  elastic  waveguides,  involving  non- 
mixed  edge  conditions  (like  stress  or  displacement),  have  remained 
unsolved.  Basically,  this  is  because  known  separation  methods  (classical 
or  integral  transforms)  do  not  "ask"  in  a  natural  way  for  the  given  edge 
information.  A  means  for  solving  some  problems  in  this  class,  focused 
on  the  semi-infinite  plate,  as  an  example,  is  presented  here.  In  the 
method  a  Laplace  transform  is  used  on  the  propagation  coordinate,  say  x. 
Exploitation  of  the  boundedness  condition  on  the  solution,  at  x-‘oo,  generates 
two  coupled  integral  equations  for  the  edge  unknowns  (displacements  and 
strains),  which  depend,  parametrically,  on  those  complex  wave  number  roots 
of  the  governing  Rayleigh- Lamb  frequency  equation  representing  unbounded 
waves.  Solution  of  these  equations  determines  the  transformed  solution  of 
the  problem,  which  can  be  inverted  through  known  techniques.  Excitation 
of  a  plate  with  a  built-in  edge  is  treated  as  an  example. 

51.  INTRODUCTION 

Attempts  at  solving  waveguide  problems  based  on  the  equations  of 
motion  from  linear  elasticity  theory,  and  involving  non-mixed  edge  conditions, 
date  back  to  Pochhammer's  classical  work  [l]^  in  1876  on  the  frequency  equa¬ 
tion  for  the  infinite  circular  cylindrical  rod.  As  Love  [2]  shows,  attempts  to 
use  this  theory  to  treat  the  free  vibration  of  the  finite  length,  homogeneous, 


Numbers  in  brackets  designate  references  at  the  end  of  the  report. 
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isotropic  rod,  with  a  separation  technique,  lead  to  a  solution  in  which  the 
normal  stress  on  the  rod  ends  vanishes,  but  not  the  shear  stress.  Only  an 
approximate  solution  could  be  written  based  on  the  argument  that  if  the  rod 
was  long  and  thin,  with  the  shear  stress  zero  on  the  cylindrical  surface,  it 
is  also  approximately  zero  on  the  ends. 

Similar  difficulties  are  evidenced  in  the  modern  work  done  in 
elastic  waveguides.  For  example,  the  integral  transform  technique  given 
by  Folk,  Fox,  Shook,  and  Curtis  [3]  for  solving  semi-infinite  elastic  rod 
problems,  involving  transient  end  load  inputs,  are  restricted  to  mixed-end 
conditions,  i.  e.  they  treat,  for  instance,  a  case  of  suddenly  applied  normal 
stress  at  the  rod  end  under  radial  (displacement)  constraint.  Simply  put, 
those  elastic  waveguide  problems,  that  have  been  solved  by  a  direct  integral 
transform  method,  have  been  geared  to  the  mixed-end  (or  edge)  conditions, 
since  these  are  what  the  transforms  "asked"  for. 

In  problems  involving  non-mixed  end  (or  edge)  conditions,  i.  e.  the 
stresses  or  displacements,  very  little  information  has  been  obtained, 
basically  because  direct  separation  methods  fail  to  yield  solutions.  A 
means  for  solving  some  problems  in  this  class,  focused  on  the  semi-infinite 
plate,  as  an  example,  is  presented  here.  In  the  method  a  Laplace  transform, 
parameter  s,  is  used  on  the  propagation  coordinate  x.  To  insure  that  a 
solution  be  bounded  at  x-*oo,  residues  of  poles  in  the  transform  of  the 
solution,  occurring  in  the  right  half  s -plane  and  corresponding  to  complex 
wave  number  roots  of  the  governing  Rayleigh- Lamb  frequency  equation,  are 
set  equal  to  zero.  This  generates  two  coupled  integral  equations  for  the  edge 
unknowns  (displacements  and  strains)  which  depend,  parametrically,  on 
these  complex  wave  number  roots.  Solution  of  these  equations  determines 
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the  transformed  solution  of  the  problem,  which  can  be  inverted  through 
2 

known  techniques  .  Excitation  of  a  plate  with  a  built-in  edge  is  treated  as 
an  example. 


f2.  GENERAL  METHOD 


Formal  Solution 


Consider  the  semi-infinite  elastic  plate  of  thickness  2h  sketched  in 

Fig.  1.  If  we  let  u,  v,  w  be  the  displacements  in  the  x,  y,  and  z  directions, 

3 

respectively,  and  assume  plane  strain  such  that  w  =  =0,  the  governing 

displacement  equations  of  motion  from  linear  elasticity  theory,  for  the 
homogeneous,  isotropic  plate,  can  be  written  as 


2  2  2  2 
Cd  uxx(x'y*t)+(cd  "  cs  )vxy*x,y,t)  +  Cs  uyy(x,y,t)  =  utt(x,y,t) 


xy' 

2  2  2  2 
c  s  vxx(x,y,t)+ (cd  _C«  K,Jx*y’t>  +  c,i  vjx-y-1)  =  vft(x,y,t) 


s  xy' 


yy 

A  V  ( 

d  yy' 


tt' 


<■) 


2  2 

where  c^  =  (X+2|J)/p  ,  and  cg  =M/p  are,  respectively,  the  dilatational  and 
equivoluminal  body  wave  speeds,  X  and  4  being  Lame's  constants.  Corre¬ 


sponding  stress- strain  relations  are 


CTx(x,y,t)  ,  ,  ( k^-  2]  , 

TW  =  ux(x- y* +  Vx’ y* l)  ’ 


<*(*.  y.O 

x+Tu 


U  {x,y,t)  + V  (x,y,t) 
A  y 


\ 


(2) 


°  (*.y.t) 

_ 

4 


=  V  (x,y,t)  +  u  (x,y,t)  , 

A  y 


and  o  =  v(a  +  a  ) 
z  x  y 


j 


The  method  developed  here  is  an  extension  of  one  developed  by  M.  Picone 
that  the  author  learned  of  through  private  communications  with  Professors 
A.  Ghizzetti  and  W.  Cross  of  the  Instituto  Nazionale  Per  Le  Applicazioni  Del 
Calcolo,  Rome,  Italy.  Picone  was  interested  in  finite  domain  problems, 
hence  used  a  finite  Laplace  transform  and  the  entirety  of  this  transform  in 
the  s-plane.  Evidently  this  work  was  not  published,  but  Picone's  related 
work  on  the  finite  Laplace  transform  [4]  was  published  in  1939. 
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where  V  is  Poisson's  ratio.  Subscripts  in  this  work,  when  associated  with 
displacement,  indicate  differentiation,  but  when  associated  with  stress 
identify  the  component  in  the  usual  way.  Initial  conditions  are  taken  as 

u(x,  y,  0)  =  ut(x,  y,  0)  =  v(x,  y,  0)  =  vt(x,  y,  0)  =  0  (3) 

and  conditions  at  x-*oo  as 


lim  \ 
x  —  co 


u,  u  ,  etc 

x 

v,  v  ,  etc 

x 


) 

r 


(4) 


Now  introducing  the  one  sided  Laplace  transforms  on  x  and  t,  the 
transforms  of  (1),  using  (3),  are 

u(s.  y,p)+  (k2-l)sv  (s,  y,p)+  (k2s2-  k  2)u(s,y,p) 

7  7  1  * 

=  k2[su(0,y,p)  +  ux(0,y,p)]  +  (k2-l)vy(0,  y,p)  =  f(y) 


vyy(s,  y,  p)+  (^-^-)sily{s,  y,  p)  + 
k 


'.z-k2 

s 


<*> 


^(s.y.p) 


=  -^•[sv(0,  y,  p)  +  vx(0,  y,  p)]  +  (^-^-)u  (0,  y,  p)  =  g(y) 


2  2  2  2  2  2 

where  k^  =  p  /c^  and  kg  =p  /cg  ,  and  the  bar  and  tilda  over  quantities 

indicate  the  time  and  space  transforms,  respectively.  It  is  important  to 

note  that  f(y)  and  g(y)  are  composed  of  the  edge  unknowns  u,  u  ,  and  v  , 

x  y 

and  v,  v  ,  and  u  ,  respectively,  these  quantities  being  the  time  transformed 
x  y 

displacements  and  strains.  Likewise,  the  transforms  of  (2)  are 

\(8,y,p)  -  /k2-2V 

+-yu-  =  su(s.y.p)  -  u(0,y,p)+  [ — j-Jv  (s,y,p) 


TT7p~ 


°v(s.y.p)  /k2.?\ 

yx  +  ^  JC8u(s,  y.  p)  -  u(0,  y,  p)]  +  vy(s,  y,  p) 


)  (6) 


«-«f*  •'wnWMah  -t 
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u(s,y,p)  =  Aa(s,p)e"ay+ A_a(s,p)eay+ Ap(s,p)e'Py 

+  A_p(8,p)ePy  +  up(8,y,p) 

v(s,y,p;  =  cpa[Aa(s,p)e'ay- A_a(s,y)eay]  +  ^[a^b, p)e'^y 

-  A_p(«.  »>)e^y]  +  vp(s,y,p) 

where 


>  (9) 


/ 


cp.  = 

OL 


Ba(s,p) 

Vs* P) 


a 

8 


Bg(8,p) 

V  Ap(.,p) 


8 

F 


(9a) 


stems  from  the  fact  that  the  two  algebraic  equations  yielding  (7a)  must  hold 
for  all  values  of  n. 

Although  certainly  not  a  restriction  to,  particular  interest  here  will 
be  in  symmetric  excitation  of  the  plate  (w.  r.  t.  the  mid-plane,  y=0).  On 
this  basis  it  is  only  necessary  to  consider  half  of  the  plate,  i.  e.  O^y  *h. 
Narrowing  further  to  plate  face-type  loadings,  as  a  special  case,  the 
boundary  conditions  at  y=h  are,  from  the  second  and  third  equations  of 
(2)  and  their  transformation  in  (6), 


^  (s.h.p)  /k2_2\r  ^  n  =*  a0F(s)G(p) 

\  V2|i  "  =  V — 2~j[su(8  »h,p)-u(0,h>p)]+Vy(8,h,p)  - \  +  =T(s,p) 

a  ( s ,  h,  p)  ^  ^  _ 

xy  — -  =  Uy(s,  h,  p)  +  sv(s ,  h,  p)  -  v(0,  h,  p)  =  0 


(10) 


where  GqF(x)G(x)  is  the  symmetric  (w.  r.  t.  the  plate  mid-plane)  normal 
stress  input  on  the  plate  face  y  =  h,  F(x)  and  G(t)  being  arbitrary,  but 
transformable,  functions,  and  0^  an  inherently  negative  magnitude  constant 
of  dimensions  force/unit  length.  At  y=0  the  conditions  are 


* 


a  (s,0,p)  _ 

-X-Y-  — - -  U  (s,  0,  p)  +  sv(s,  0,  p)  -  v(0,  0,  p)  =0 

v(s,  0,  p)  =0  J 

Substitution  of  (9)  in  (10)  and  (11),  determines  the  coefficients  A 


±a,  ±p ' 


They  are 


Aa=  A-a=  A(8,p)  = 


A^P) 

R(s,p) 


BN(8»  I  ) 

Ap=A.p  =  B(.,P)=-Ni7?r 


where 


2A^(s,p)  =  -  s£k2(2s2-kg2)sinhph  •  I+2spcoshph  • 

2BN(s,p)  =  -p[2k^sa  sinhah  •  I  -  (2s2-ks2)cosh ah  •  J~! 

/  2  2  \2  2 
R(s,p)  "  \2s  -kg  J  coshah  si.ihph+ 4s  ap  sinhah  cosh  ph 


k  o  k 

8 


,  I  (2s2-k  2) 

I  =  T(s,p)+ — j]  \-jr  - - — —  sinha(h-y ')  +  2p  sinhp(h-y  ')  f(y#) 


+  £(282-ks2)cosha(h-y ')  -  2s2  cosh P(h-yr  )Jg(y  ^dy ' 
+  (— -7^)u  (0,  h,  p) 

J  =  — -j  J  /^s2  cosha(h-y  ')  -  (2s2-kg2)cosh P(h-y 7 )^|f (y  ') 


ks  0 


s 

+  ^p-[2ap  sinha(h-y ')  +  (2s2-k82)sinh P(h-y  ' )]g(y  'jjdy ' 


W0,  h,  p) 
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where  it  should  be  pointed  out  that  and  contain  the  edge  unknowns, 
through  f(y)  and  g(y),  and  the  corner  displacements  u(0,h,p)  and  v(0,h,p), 
in  I  and  J.  Formally,  however,  the  transformed  solution  on  t  can  be 
written  as 


and 


u(x,y,  p)  [uc(s,y,p)  +  up(s,y,p)]e8xds 


v(x,  y.p)  '  2iri  Jg  [vc<8*  Y»  P) +  vp(a.  y,  p)]e8Xds 


\ 


J 


(13) 


where  the  complementary  functions  are  given  by 


and 


uc(8»y»P)  =  2[A(s,  p)coshCty  4- B(s.  p)coshpyJ 


vc(s,y,p)  =  2[|A(s,p)sinhay-^-B(s,p)sinhpy] 


where  Brg  is  the  Bromwich  contour  in  the  right  half  s -plane. 


Exact  Inversion  Procedure  for  Steady  Propagation  Case 

Since  the  integrands  in  (13)  are  even  functions  of  both  a  and  p,  they 
are  singled  valued  functions  of  s  everywhere  in  the  s-plane,  for  an  arbitrary, 
but  fixed,  value  of  p.  This  reflects  the  finiteness  of  the  domain  in  y.  Setting 
p=iu),  (13)  represents  the  formal  solution  for  the  steady  wave  propagation 
case  in  the  present  problem.  Further  if  s=ix,  R  becomes,  in  this  case, 

R(s,  p)  |  =ijjx^-p/^)  cosa,hsinp/h+ 4x^a,p/ sina,hcosp,hJ  (14) 

p=iu) 

8  =iH 

where 
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U)  and  K  being  frequency  and  wave  number,  respectively.  The  bracket  on 
the  right  hand  side,  set  equal  to  zero,  may  be  recognized  as  the  classical 
Rayleigh- Lamb  frequency  equation  for  symmetric  straight-crested  har¬ 
monic  waves  in  an  infinite  elastic  plate,  where  a  and  p'  are  the  thickness 
wave  numbers  in  the  plate  for  dilatational  and  equivoluminal  waves, 
respectively. 

Clearly,  then,  the  roots  s^(tu)  =  iK^(tu)  of  R  in  (14),  which  are  known 
to  be  simple  zeros,  infinite  in  number,  are  the  modes,  or  branches,  of 
wave  propagation  in  the  present  problem,  at  least  for  the  steady  case. 
Through  the  efforts  of  Holden,  Mrndlin,  Onoe,  and  others,  the  nature  of 
these  modes  are  now  well-known.  A  recent  study  by  Onoe,  McNiven,  and 
Mindlin  [5]  on  the  closely  related  Pochhammer  frequency  equation^  for 
axially  symmetric  waves  in  an  infinite  circular  cylindrical  rod  gives  a 
detailed  account  of  the  main  features  of  these  spectra.  As  Fig.  6  of  [5] 

(Fig.  2  here)  exhibits,  it  is  known  that  for  real  frequency  (dimensionless 
Q  in  Fig.  2)  there  are  an  infinite  number  of  branches  (or  roots),  each  of  which 
is  a  continuous  one  to  one  relation  between  frequency  and  wave  number 
[dimensionless  C=5+iTl  in  Fig.  2,  which  is  essentially  Hh  in  (14)]  over  the 
frequency  domain  from  oo  to  -oo.  Over  this  domain  the  wave  number  may 
be  real  or  take  on  imaginary  and  complex  values,  as  Fig.  2  exhibits.  It  is 
easily  shown  that  £(0)  and  -£”(£})  are  conjugate  roots  of  (14),  the  latter  being 
a  reflection  in  the  plane  ?=0.  These  roots  are  indicated  in  the  Figure.  Now, 

since  £  occurs  only  as  Q?'  in  (14),  -C(Q)  and  £(Q)  are  also  roots,  these  being 

- - - ■ — - — — - - - 

J  Equations  (la)  and  (lb)  in  [6]  exhibit  the  closeness  in  form  of  the  two 
equations. 
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reflections  of  the  former  two  in  the  plane  11  =  0.  So  for  ft^O,  a  complex 
segment  of  a  branch  has  images  in  both  the  planes  §  =  0  and  T)=0,  an  imagin¬ 
ary  segment  an  image  in  the  plane  T|=0,  and  a  real  segment  an  image  in  the 
plane  5=0.  Fig.  2  offers  a  scheme  for  following  a  particular  branch 

(numbered)  continuously  through  its  various  segments  as  ft  progresses 

4 

from  oo  to  zero  or  vice  versa  .  It  follows  from  the  foregoing  discussion 
that  once  the  edge  unknowns  are  determined,  inversion  of  the  comple¬ 
mentary  functions  in  (13)  can  be  accomplished  with  simple  residue  theory. 
Assuming  then,  the  inverse  of  u^  and  v^  can  be  found,  this  would  complete 
the  solution  for  the  case  of  steady  propagation  in  the  present  problem. 


Exact  Inversion  Procedure  for  Transient  Case. 

Extension  to  the  transient  solution  requires  inverting  the  integrals 


in 


u(x,  y,  t) 


u(x,  y, p)eptdp 
r 

P 


v(x, y,t) 


1 

2tt  i 


4 


v(x ,  y, p)eptdp 
r 

P 


> 

/ 


(15) 


where  u  and  v  are  given  by  (13),  Br^  being  the  Bromwich  contour  in  the  right 
half  p-plane.  If  now  it  is  assumed  the  roots  s^(uu)=iK^(uu)  are  analytically 
continuable  from  y  on  the  imaginary  axis,  i.  e.  roots  of  the  right  hand  side 
of  (14),  to  the  corresponding  set  s^(p)  =  iH^(p)  on  the  Bromwich  contour  Br^, 
(15)  can  be  inverted  by  termwise  contour  integration  of  the  terms  in  the 
series  resulting  from  the  inversion  of  (13).  This  method  of  inverting  the 


4 - 2 - 

Since  ft  occurs  as  ft  in  (14),  -ft  will  also  yield  the  spectra  shown  in  Fig.  2, 
i.  e.  the  reflection  in  the  plane  ft=0.  Such  modes  usually  play  a  part  in 
deriving  the  solution  for  a  problem  (like  the  present  one).  The  reflection 
principle,  however,  reduces  the  solution  to  dependence  on  ft  s  0  only. 
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double  integral  transform,  i.  e.  inversion  of  spatial  transform  first,  was 
discussed  further  and  used  by  Lloyd  and  Miklowitz  [7]  in  related  work. 

It  is  based  also  on  the  assumption  that  there  are  no  singularities  in  u  and 
v  in  the  right  half  p-plane,  Rep  >  0.  Clearly  such  singularities  would  be 
ihconsistent  with  the  fact  that  bounded  long  time,  or  static  solutions,  for 
the  present  problems  do  exist.  It  follows  that  the  contour  integration  can 
be  carried  out  along  the  imaginary  axis  in  the  p-plane,  yielding  a  transient 
solution  in  the  form  of  a  series  of  integrals,  each  term  of  which  represents 
a  branch  of  the  frequency  equation  depicted  in  Fig.  2. 

Boundedness  Condition;  Integral  Equations  for  Edge  Unknowns. 

It  follows  from  the  discussion  of  the  roots  °f  R  in  (14)  and  their 

continuation  to  Kjc(p).  and  Fig.  2  depicting  them,  that  there  are  two  sets  of 

complex  branches,  i.  e.  s.(p)  =  iH  .(p),  and  their  conjugates  s  (p)  =  iH.(p), 

J  J  J  J 

that  satisfy 

Re  sR(p)  >  0  ,  (k  =  j)  (16) 

which  in  terms  of  Fig.  2  is 

Re  [i(5  +  iri)]  >0 
or 

T|  =  Im  C  <  0,  (16a) 

i.  e.  reflections  of  the  complex  branches  in  Fig.  2  in  the  plane  r)=0.  It  is 
clear  that  these  roots  of  R,  which  lie  in  the  right  half  s-plane,  would  give, 
through  the  residue  evaluation  of  (13),  exponentially  unbounded  waves  at 
x-^oo,  hence  violating  the  condition  (4).  It  follows  that  these  unbounded 
waves  can  be  eliminated  by  requiring  the  residues  in  (13),  associated  with 
them,  be  zero,  i.  e. 
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AN(s,p) 


=  BN(s,p) 


=  0 


(17) 


S=Sj(p) 


8=8j(P) 


where  s.(p)  are  the  two  sets  of  roots  of  R(s,p)  in  (12)  satisfying  (16). 

J 

Substitution  of  and  from  (12)  into  (17)  gives  two  equations  for  I  and 

(18) 


J,  for  each  sj(p)*  Compatibility  of  these  equations  requires  that 


J.  -  P.I.  =  0 
J  J  J 

where  J.  and  I.  are  J  and  I  in  (12)  for  s  =  s.(p),  and 
J  J  J 


P. 

J 


2k2s  .a.  sinha.h  k2(2s  2-k  2J  sinhfl.h 
-  J  J  J  _  ^  J  8  /  Kj 

~(2s2-k2)—  ^  ’  2s.p.  coshji.h 


cosha.h 
J 


a.  =  /k2-s2'  ,  p.  =  7k2  -  s?' 

j  V  d  j  V  s  j 

Now  arguing  on  the  basis  of  Lerch's  theorem^,  that  for  the  present  p  can 
be  assumed  to  be  real,  and  noting  again  the  conjugate  nature  of  the  roots 
Sj(p)  satisfying  (16),  the  boundedness  condition  (18)  can  be  expanded  into 


c  r 

Re 

.  rh  1 

/ 

Im 

ks  0 

l  J 

>  L 

28.  cosha.(h-y  ) 
J  J 


s.(2s2-k2) 

Vi  Jv  j 


2 - P-  sinha  (h-y  ) 

k  a.  J  J 
J 

2s. p. 


-  ^28j2 -kg2)  cosh (3^ (h-y*)  -  ^  J  Pj  sinh |3.(h-y ') 

2k2s.a  sinha.(h-y,)-(2s.2-  k  2)P.  cosha.fh-y*) 

J  J  J  7  '  J  8  /  j  j' 

k 2  8 .  (2  s  2  -  k  2) 

+__Aj — lL 


Hy) 


sinhp.(h-y')4  2s^  P.  cosh Pj(h-y')| g(y') 


-  v(0,h,p)-P. 


.  2  ' 

T(s;,  P)  +|  — j — j 
k  1 


Jl  '  J 


J 


u(0,  h,  p)  )  =0 


dy 

(19) 


An  inverse  Laplace  transform  is  unique  provided  the  transform  is  known 
for  all  real  values  of  p  for  p>p^  (see  Widder  [8]  ). 
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(19)  are  two  coupled  integral  equations  for  the  edge  unknowns  for  each  set 

of  parameters  p  and  s.(p).  For  a  particular  problem  the  edge  conditions 

J 

would  reduce  the  number  of  edge  unknowns  in  (19),  and  assuming  a  solution 
of  these  equations  can  be  found  for  the  remaining  edge  unknowns,  this 
would  complete  the  formal  solutions  given  by  (13)  and  (15)  as  will  soon  be 
demonstrated. 

§3.  PROBLEM  OF  PLATE  WITH  BUILT-IN  EDGE 

As  a  first  example  the  case  of  a  plate  with  a  built-in  edge  was 
treated.  As  illustrated  in  Fig.  3  the  plate  is  excited  symmetrically  by  two 
suddenly  applied  normal  line  loads  on  its  faces,  a  distance  a  from  the  edge. 
This  particular  loading  is  expressed  in  the  first  of  (10)  through 

o  (x.h.t)  a  6(x-a)G(t) 

Utu  <20a: 

where  ^(x-a)  is  the  symmetrical  delta  function,  so  that  T(s,p)  in  (10),  and 
later  equations,  is 

a  G(p)e‘sa 

T(s,p)  =  (U2U) 

The  edge  conditions  are  given  b/ 

u(0,  v,  t)  =  v(0,  y,  t)  =  0 
u  (0.  y,  t)  =  v  (0,  y.t)  =  0 

so  that 

u{0,  y,  p)  =  v(0,  y,  p)  =  0 
u  (0,  y,  p)  =  v  (0,  y,  p)  =  0 

(21),  it  may  be  noted,  reduces  the  unknowns  in  the  two  equations  (19)  to  two, 
u  (O.y.p)  and  vx(0,y,p),  hence,  apparently,  a  well-posed  and  simpler 


(20b) 


(21a) 


(21b) 
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problem  in  the  present  analysis  than  the  free  edge  (or  edge  loaded)  prob¬ 
lem.  In  the  latter  one  would  have  to  deal  with  tht  four  unknowns  u(0,y,p), 
ux(0,y,p),  v(0,y,p),  and  TjO,  y,  p). 

Determination  of  Edge  Unkowns  for  Long-Time  Solution. 

During  the  course  of  this  work  the  author  learned  of  a  paper  by 

Benthem  [9]^  who  exploited  a  similar  boundedness  condition,  to  that  used 

here,  to  solve  related  static  problems  in  strips  governed  by  the  biharmonic 

7 

equation  on  the  Airy  stress  function  for  plane  stress  .  Of  particular  interest 
in  Benthem's  technique  is  his  use  of  eigenfunction  expansions  (Fourier  series 
in  his  problems),  and  singular  terms  (where  needed),  to  represent  the  edge 
unknowns.  This  scheme  reduces  the  boundedness  condition  (in  his  work)  to 
a  finite  set  of  algebraic  equations  for  the  coefficients  of  the  Fourier  series 
and  the  singular  terms.  With  numerical  evolutions,  Benthem  shows  that 
only  a  very  limited  number  of  the  coefficients  need  be  determined.  Benthem's 
comparisons  of  his  results  with  those  in  the  early  work  of  Knein  [ll]  for 
the  problem  of  the  strip,  or  plate,  with  a  built-in  edge,  which  the  latter 
obtained  by  other  means,  show  very  good  agreement.  It  was  therefore  of 
interest  to  see  if  Benthem's  technique  could  be  extended  to  finding  solutions 
of  (19),  and  in  particular,  for  the  problem  of  the  built-in  edge  being  treated 
here. 

Since  on*  would  expect  the  static  solution  as  the  long-time  limit  in 

the  present  problem,  it  seemed  reasonable  to  formulate  representations  for 

^  The  author  is  indebted  to  his  colleague  Professor  James  Knowles  for 
pointing  out  Benthem's  work. 

y 

Benthem  points  out  that  Doetsch  [10],  in  a  strip  problem  governed  by 
Laplace's  equation,  employed  a  boundedness  condition  in  a  manner 
related  to  that  in  Benthem's,  and  hence  the  present  work. 
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the  edge  unknowns,  with  expansions  similar  in  form  to  those  used  by  Benthem 
for  his  problem  involving  a  built-in  edge,  and  derive  the  long-time  solution. 
On  this  basis,  it  was  assumed  that 


f(y)  =  k  ux(0,  y,  p)  = 


-k‘ 


a0<P> 

K‘-£)q 


l  a„<p> 

n=l,3,5,.. 


niry 

cosTh 


(22a) 


M  /1'i-m  ,  Ca0<P>r  1  /.  y\l  ,  V  bn(p)  •  niry 

sM =  ^vx<°*  y>  p) =  nyq  •  V  ■  t)  +  l  -pr  8in~z£ 

V  n  /  n=2,4,6,.. 


(22o) 


These  are  similar  in  form  to  Benthem' s  equations  (6.  1)  for  the  edge  normal 
and  shear  stress.  In  comparing  the  two  it  may  be  noted  whereas  he  had 
constant  coefficients,  here  a^,  a^,  and  b^  are  functions  of  the  parameter  p. 
The  power  q  of  the  corner  singularities  in  (22)  is  taken  to  be  that  found  in 
the  eigenfuncti^^  problem  for  the  clamped-free  wedge.  This  was  treated  by 
Knein  [  1 1  ]  for  the  right  angle  wedge  or  corner  (case  in  [9]  and  here),  and 
by  Williams  [l2]  for  more  general  angle  corners.  As  these  works  show  q 
depends  only  on  Poisson's  ratio  for  a  fixed  corner  angle.  Hence,  although 
Knein  treated  a  case  of  the  built-in  edge  (and  plane  strain)  for  a  uniform 

edge  compression,  whereas  here  lateral  loads  are  involved,  a  particular 

8  |  ^xv I 

q  holds  for  both  problems  .  Knein  gives  curves  for  q,  and  <■  along 

'  ax  ' 

x=0  (near  the  corner)  which  determines  C  in  (22b),  as  a  function  of 
Poisson's  ratio  v.  Following  Knein,  v  is  taken  to  be  0.  2433  here,  which 
makes  q  very  close  to  1/4  ,  and  C  =  0.  740. 

Then  by  substituting  (22)  in  (19),  and  with  the  aid  of  the  integrations 


g 

The  same  q  is  also  applicable  in  Benthem's  problem,  which  involves  an 
applied  uniform  tension  load  over  y  at  x-*oo.  Since  Benthem's  case  is 
one  of  plane  stress,  Poisson's  ratio  must  ^e  different. 


-16- 


/jit} 


n-l 

(-)^  a.  s'.nha.h 
'  '  2h  j  j 


N 


n 


a.  cosh  a. h 
J  J 


>,  11=1,3,5,. .  (23a) 


and 


hr 

sinh 


cosh 


{^(h-y^sin^Xdy'  =  -U 


a: 

j 


^  sinh  a.h 

coshah] 


>.  n=2,4,6,.  .  (23b) 


and 


(l-y'/h){  )  {a  (h-y')j  dy'  =  -j-J 

sinh  I  J  a^  h  ] 


a.h  sinh  a.h  -  cosha.h  +  1 
J  J  J 

a.h  cosha.h  -  sinha  .h 
J  J  J 


and 

ib  h--')  * 


■ 


a. 

J 


'a. 

J 


(24a) 


(24b) 


n  2  2  2  2 

where  A^n=  a^  +  n  w  /4h  ,  and  the  aid  of  the  duplicate  set  of  integration 

to  (23)  and  (24),  where  (3.  and  B.°=  p.^+ n^ir^/4h^  replace  a.  and  A.n, 

J  J  J  J  J 

respectively,  (10)  can  be  reduced  to  the  algebraic  equations 

n-l 


n=l,3,5,.. 


r  \ 

Re 

r 

>  j 

Im 

n 

,7  n. 


+  )  (-)  N.  (s.,  p)b  (p)  -  Q.(s 

J  J  11  J 

o=2,4,6,. . 


i(8j  'P)j 


=  0 


(25) 


where 
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(25)  can  be  solved  for  the  unknowns  a^p),  an(p)»  and  bn(p).  That 
is,  for  a  certain  number  of  unknowns,  a^p),  aj(p),  a^(p),-",  b^(p), 
b^(p),  •  •  •  ,  a  matching  number  of  s^(p)  (which  are  infinite  in  number; 
see  Fig.  2)  are  available  to  give  a  sufficient  number  of  equations  from 
(25)  to  solve  for  these  unknowns.  Benthem  shows  in  the  static  case  by 
progressively  taking  more  and  more  unknowns  (constants  like  a^'s  and 
b^'s  here),  he  gets  convergence,  numerically,  to  an  a^  and  b^  value, 
where  N  is  a  relatively  small  number.  As  he  points  out  this  would  have 
to  be  the  case  since  the  series  representations  in  (22)  ,  obviously,  are 
not  being  called  upon  to  represent  the  singular  terms  there,  i.  e.  if  they 
were,  then  one  would  expect  the  convergence  to  be  quite  slow.  It  follows 
that  in  the  long-time  solution  for  the  present  problem  analogous  behavior 
could  be  expected. 


-18- 


The  long-time  solution  can  be  derived  from  (15)  by  making  use  of 
the  underlying  Tauberian  theorem,  i.  e. 


lim  ( 
t-oo 

r  ] 

u(v,y,t) 

v(x,y,t) 

>  =  lim/ 

p-*0 

1 

2ni 

r 

Br 

P 

e-  _ 

u(x,  y,p) 

v(x,  y,p) 

“  •- 

N 

ept dpi 

„  J 

J 

(26) 


Hence,  the  solution  of  (25)  must  be  obtained  for  p-*0  in  order  to  determine 
u  and  v  in  (26).  The  roots  s^(p)  of  R(s,p)  in  (12)  can  be  found  by  evaluating 


=  lim  (  c 

p-0 


2  3Rl 


s 


(27) 


3p 


which  gives 


.  2 
1 8 

^(sin  2sh+  2sh) 


Hence  the  s.(p)  are  selected  from  the  zeros  of 
J 


(28) 


r(s)  =  sin  z  +  z  ,  z  =  2sh  =  x'  +  iy' 


(29) 


satisfying  (16),  i.  e.  Res^X).  The  zeros  of  r(s)  are  an  ordered  infinite 
set,  corresponding  to  the  piercing  points  of  R  in  the  plane  £1=0  as  Fig.  2 
indicates.  r(s)  is  a  well  known  function  in  the  analyses  of  elastostatic 
plate  problems.  Hence  its  occurrence  here,  and  in  Benthem's  work  is  not 
surprising.  Robbins  and  Smith  [13]  give  the  first  ten  roots  to  6  decimals, 
and  the  asymptotic  behavior  of  large  roots.  The  behavior  of  the  slopes  of 
Sj(p)  can  be  found  by  noting  that  along  these  branches 

dR  =  "5pdp+  T7d8  =  0 


and  therefore 


lim 


p-°  P-57 


(30) 


-19- 


which  from  (28)  is 

«3i> 

p-0  K  p  r  (s) 

since  r(s)  =  0,  and  r  (t  )  £  0  ,  the  prime  denoting  differentiation  w.  r.  t.  s. 

It  follows  that  ds/dp  in  (31)  must  vanish  in  the  limit,  i.  e.  these  branches 
are  normal  to  the  plane  0=0.  It  follows  that  the  zeros  of  r(s)  are  a  good 
approximation  to  the  s^(p)  for  a  range  of  p,  small  but  >0.  This  is  impor¬ 
tant  to  the  validity  of  (26),  hence  the  long-time  solution  being  sought.  It 
should  also  be  pointed  out  that  (28)  directly  proves  the  existence  of  the 
zeros  of  R(s,p),  at  least  over  the  domain  of  small  complex  p,  that  was 
argued  earlier  on  the  basis  of  the  analytical  continuation  of  these  zeros, 
generally,  from  the  domain  of  imaginary  p  [Re  discussion  after  (15)1. 

With  the  Sj(p)  determined  from  the  zeros  of  r(s)  it  remains  to 
approximate  the  coefficients  in  (25)  for  p  small,  and  then  solve  for  the 
unknowns.  Expanding  a.,  p^,  and  quantities  dependent  on  these,  for  small 
p,  there  results 


a.(s  ,p)=is  1 - ?  +  0(p  )  ,  p(s  ,p)=is  1 - 5  +  0(p  ) 


ih  l  a  ^ 

=  sinhis.h--^ — <  /cosh is .h  +  0(p  ) 

J  ^s;  y  J 


cosh 


=  cosh  iSjh  -  )sinhis.h+ O(p^) 


which  can  be  used  to  obtain 
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where  S^g  and  C^g  are  given  by  the  second  integral  in  (24b),  where 

J  J 

a.  =  Q. = is. ,  and 
J  J  J 

Y.  =  [  C^cosh  is  .hC  dC 
j  o  J 


6.  =  [  C^sinh  is.hC  dC 

1 S  .  J  j 

j  o 
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p-»0 


Mjn(“j’P»  - 
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2hUn 
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+  0(p4) 


lim  N.n(s.,  p)  = 

p-*0  J  J 


inirx. 
is . 

- 

2hk  Un 
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2 

s 


+  0(p4) 


\  (32) 


and 


•s.a 


lim  Q  (s  ,p)  =  i(a  /n)e  j  x  G(p)k^  [  1+0  (p2)] 
P"*0  J  J  j 


y 


where 


x. 

IS  . 

J 


sinhis.h 

,  J.  • 
cosh  is  .h  ’ 
J 


2  2 

,  TT  n  n  it 
and  U.  =  - 

J  41^ 


Since  the  real  and  imaginary  parts  of  the  leading  terms  in  (32) 
behave  as 


their  substitution  in  (25),  shows  that  a  solution  of  the  latter,  say  by 
Cramer's  rule,  will  have  the  behavior 


limfa0(p),  aj(p),  a^(p),  .  .  .  b^p),  b4(p),b6(p), .  .  .]~G(p)  (33) 

p-*0  L 

(33)  is  an  important  result  since  it  shows  the  long-time  behavior  of  the 
unknown  strains  at  the  edge  (x=0),  ux(0,y,t)  and  vx(0,y,t),  according  to 
(22),  is  the  same  as  the  input  loads  on  the  plate  faces.  Further  if  G(t)  is 
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set  equal  to  the  step  H(t),  it  is  clear  these  edge  unknowns  have  the  same 
form  as  in  Benthem's  work,  i.  e.  the  unknowns  become  constants.  (33), 
and  (32),  show  the  numerical  problem  remaining  to  completely  define  the 
unknowns  is  simple.  One  would  evaluate  forms  such  as 

a3(p) .  bnl 


where  the  R.  H.  S.  depends  only  on  n  and  s^,  and  not  on  p  [ s ^  are  the  roots 
of  r(s)  in  (29)1,  and  use  Benthem's  numerical  convergence  procedure  for 
the  unknowns,  i.  e.  repetition  of  the  solution,  with  a  larger  number  of 
unknowns,  until  convergence  to  the  final  value  of  each  unknown  is 
reached.  Benthem's  example  shows  convergence  is  reached  quickly. 

Except  for  carrying  out  this  numerical  work,  the  analysis  given  here  makes 
it  possible  to  define  the  formal  long-time  solution  (26).  In  addition  to  the 
very  long-time  (static)  solution,  (26)  will  yield  the  low  frequency-long 
wave  dynamic  contribution  from  the  lowest  real  mode  (see  Fig.  2)  that 
occurs  earlier.  Since  p  is  also  small  for  this  contribution,  it  can  be 
derived  from  (26)  because  it  is  compatible  with,  and  hence  can  be  based 
on,  the  approximations  contained  in  (32)  and  (33).  The  numerical  work 
required  to  define  the  unknowns,  and  the  study  of  the  low-frequency  and 
static  contributions  to  the  long-time  solution  (26),  are  planned,  and  will 
be  reported  in  a  sequel  to  this  work. 

Formal  Long-Time  Solution. 


With  the  transformed  edge  unknowns  u  (0,y,p)  and  v  (0,y,p) 

X  X 

determined  for  small  p,  through  the  analysis  given  in  the  previous  section, 
the  formal  solution  for  long-time  given  by  (26)  can  be  determined,  i.  e. 


w 


* 


r- 


V 
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u(x,  y,  p)  and  v(x,y,p)  in  (26)  can  be  determined.  Returning  to  (13),  the 
integrands  there,  corresponding  to  (22),  can  be  obtained  by  substituting 
the  latter  into  (8)  and  (12)  and  performing  the  integrations  indicated.  The 

-~w  — 

resulting  double  integral  transforms  for  u  and  v c  are  given  by 
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Sa  and  Ca  being  the  integrals  formed  from  (24b)  by  dropping  the  sub  j 
(Sp  and  Cp  obtained  from  these  by  replacing  a  with  P). 

The  double  integral  transforms  for  u^  and  are 
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where  in  getting  u^  and  v^j  ,  use  has  again  been  made  of  (23),  replacing 

h  with  y  and  a.  with  a,  etc.  It  should  be  noted  that  in  (35)  and  (36)  the 
J 

subscripts  1,  and  2,  correspond  to  the  trigonometric  and  singular  terms 

in  (22),  respectively,  and  the  subscript  3  to  the  linear  term  in  (22b). 

These  subscripts  help  to  identify  related  terms  in  (35)  and  (36),  i.  e.  in 

u  and  u  and  v  and  v  .  (36)  shows  that  all  u  .  and  v  .  are  zero  at  y  =  0  in 

c  p  c  p  '  pi  pi  7 

accord  with  (8)  from  which  they  stem.  It  may  be  noted  that  the  super  zero 

=»  0  =*  0 

terms  ucj  and  u^  are  equivalent,  except  fcr  sign,  to  the  coshay  and  coshpy 

terms  in  u^j  and  respectively,  and  hence  their  sum,  in  each  case, 

which  pertains  to  y>0,  does  not  contribute  to  u.  Indeed  the  super  zero 

means  that  these  terms  are  in  the  solution  for  u  only  at  y=0.  In  the  case 

of  u^  no  super  zero  term  appears.  Its  integral  form  precludes  separating 

9 

such  a  term.  A  term  of  this  nature  is  involved,  however  .  The  super  R 

terms  are  so  labled  because  of  the  1/R  in  them,  and  they  contribute  to  the 

solution  for  all  y.  In  the  case  of  v  .  and  v  .  like  cancellation  is  involved 
7  ci  pi 

for  y>0,  but  since  v  =  0  at  y=0,  there  are  no  super  zero  type  terms,  and  except 

for  v  and  v  ,,  all  terms  have  been  reduced  (to  the  forms  similar  to  u  . 
c2  p2  '  ci 

and  u  .  for  y  >  0). 

pi 

Important  is  the  fact  that  (35)  and  (36)  have  resulted  only  from  the 
assumed  form  of  the  unknowns  given  by  (22),  i.  e.  they  are  not  directly 
dependent  on  the  assumption  of  small  p.  For  use  in  the  long-time  solution, 
however,  they  will  have  to  be  approximated,  for  small  p,  to  be  consistent 
with  the  approximations  given  by  (32)  and  (33)  for  the  unknown  coefficients 
3q(p),  an(p)>  and  bn(P)>  they  contain.  The  implication  here  is  that  (22)  may 


ouch  a  term  is  recovered  in  the  near  field  solution  derived  later,  and  shown 
to  exist  for  only  at  y  =  0. 
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be  valid  for  the  propagation  of  other  waves,  e.  g.  harmonic  waves  for 
certain  single  frequencies,  calculated  from  (13).  This  will  be  investigated 
in  future  work. 


Approximation  for  the  Near  Field. 


The  near  field  solution  (x<<  a)  can  be  written  by  applying  the  under¬ 
lying  Tauberian  theorem  for  this  case  to  (13),  with  (35)  and  (36),  expanding 
terms  in  the  latter  for  |  s h |  > >  1 ,  with  |k  / s  |  <  1,  making 


=  -is  1 


k2/2s2 

k2/2s2 

s 


4 

*°  ^ 


It  should  be  noted  that  since  large  |  s|  here  is  necessarily  on  Brg,  (37) 
shows,  assuming  again  that  p  is  real,  that  a  and  (3  are  real  and  large 
through  the  value  of  Im  s.  Note  that  for  very  large  |  s  |  ,  a=p=Im  s 

In  accord  with  (35)  and  (36)  the  near  field  solution  for  y=0,  and 
y  >  0  differ  [see  remarks  after  (36)1.  Because  of  the  singular  terms  in 
(22),  y=h  must  also  be  given  special  consideration.  Focusing  first  on  the 
"interior"  solution,  0<y<h,  the  contributions  of  uc  j  and  j  to  the 
transformed  long-time,  near  field  solution  are 


=  o[e1 


A3]  I 

— [  lj  3|  5|  •  • 

+  I  <-»r2Kb„<P» 

n=2,4,6,.. 


(38a) 


Since  Im  s  has  a  change  in  sigjn  from  lower  to  upper  half  of  Brs  ,  so  do  a 
and  p.  However,  since  u  and  v  are  even  in  a  and  p,  it  follows  the  large  is  hi 
approximation  holds  over  all  of  Brg  in  the  usual  way.  In  the  work  that  fol¬ 
lows  s  is  chosen  on  the  upper  half  of  Brg  . 
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which  is  based  on  the  fact  that  R=-(kd2-kg2  )s2e  2lsh/2  for  large  |sh|. 
corresponding  terms  u^j  and  v^  give 


The 
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n=l,  3, 5, .  . 


pl 
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Now  since  S  and  C  may  be  written  in  the  form 

a  a  ' 


(sj 

rh 

ea(h-y')  ?e-a(h-y')’ 

raJ 

2h(l-y'/h)1/4 
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for  large  a  these  become 

S_ 
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e  |  e  7  dy 
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(38b) 


(39) 


(40) 


Since  a  is  a  real,  positive,  large  parameter  in  (40),  according  to  (37),  this 
integral  may  be  approximated  by  Watson's  Lemma  [  1 4 D ,  with  the  result 


(41) 
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with  an  equivalent  set  for  and  C^.  (41)  may  now  be  used  to  approximate 

u  and  v  0  with  the  results 
c2  c2 
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A  parts  integration  shows  that  the  corresponding  ♦'.rms  up2>  and 
the  same  order  of  approximation,  are 


Expanding  and  vc^  it  is  found  that 
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5^^  and  v  ^  are  given  by 
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(43b) 


^  ~  S  d 

Since  and  v  ^  contain  the  shift  operator  e  they  contribute  nothing  to 
the  solution  in  x<a,  and  hence  nothing  to  the  near  field  solution. 

From  (38)  it  follows  that 


S1  =  "~Z  1  an(p)COSThi+0('T) 


n=l,3,5,.. 
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and  from  (42) 


11 


and  from  (43) 


1 1 


an<P) 
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(45) 


(46) 


It  follows  from  (44)  to  (46)  that  the  time  transformed  long  time-near  field 
solution,  valid  for  the  plate  interior,  0<y<h,  is  given  by 


u(x,y,p)~  - 


v(x,  y,  p) 


an(P)  r1 

.  '  .  V.1M+  l  an(p)COS^lf 

LMl-y/h)  n=i,  3, 5, . . 
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Ca0(P)f  i 


|(  i -y/h) 


T7?'(1_h  )/+I  bn<P>8iniS 


rt=2,4,6,.. 


/  (47) 


The  corresponding  solution  for  y=0,  is  written  by  first  noting  that 

a  ^  ai  ^  o 

u  ,  v  ,  and  v  are  all  zero,  and  hence  both  u  .  and  u  .  must  be  evaluated 
p  p  c  ci  ci 

— w.  — w.  — w.  Q 

in  (35)  to  calculate  u(s,  0,  p)  =uc(s,  0,  p).  It  is  found  by  approximating  u^ 

in  (35),  and  adding  it  to  u  ^  in  (38a),  evaluated  at  y=0,  that 

*  where  cancellation  of  e  isy  terms  occurred  to  the  order  indicated.  From 
the  pattern  found  in  these,  it  is  reasonable  to  expect  that  the  same  pattern 
of  cancellation  would  result  to  any  arbitrary  higher  order. 

When  v3  is  used  in  connection  with  the  y=h  solution,  this  term  is  replaced 
I  s3). 


by  0(  1  /  s-5) 


12 
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Uj(s,0,p)  =u^j(s,0,p)  -  — j  ^  an(p)  +  o(  — j  ) 


(48) 


n=l,3,5,.. 


— w 

Recalling  that  u^2  has  both  super  R  and  zero  contributions  in  it,  the  first 
of  (42a),  evaluated  at  y=0,  gives 


u2(s, 0,  p)  =  uc2(s,  0, p)  = 


(49) 


13 


0  . 


And  by  approximating  u^  in  (35),  and  adding  it  to  u^  in  (43a),  evaluated 
at  y=0,  there  results 


CaQ(p) 


u3(s,0,p)  =uc3(s,0,p)  - - £ 


1  - ) 
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;rr'2 

2k  s 
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(50) 


Adding  (48),  (49),  and  (50),  and  again  noting  that  v  =  0,  the  time  transformed, 
long  time-near  field  solution  for  the  plate  mid-plane,  y-0,  is 


u(x,  y,  p) 


v(x,  y,  p)  =  0 


a0(p)  r 

-Tr-+  l  an<p) 

ml, 3, 5,..  _ 


\ 


(51) 


) 


It  may  be  noted  that  (51)  can  be  obtained  from  (47)  by  letting  y-*0.  One 
would  expect  this  from  continuity  arguments.  However,  the  process  used 
here  to  derive  (47)  and  (51)  is  instructive  in  showing  that  (47)  is  contributed 
by  the  particular  integrals  of  u  and  v,  whereas  the  firs'  of  (51)  comes  from 
the  complementary  function. 

The  near  field  solution  is  not  singular  when  y=h.  This  is  clear  from 
(35)  and  (36),  if  it  is  noted  that  u^2  and  v  2  involve  integrable  singularities 
only.  This  means  that  the  singular  terms  in  (47)  are  not  present  in  the 

The  additional  1/2  in  these  terms  stems  from  e1S^  in  coshay  which 
survives  for  y  =  0  only. 
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solution  for  y=h.  Now  since  all  other  terms  in  (47)  vanish  at  y=h  the 

14 

near  field  solution  now  comes  from  the  order  terms  in  (44)  and  (46)  , 


i.  e. 


u(s,  h,  p) 


M4) 


(v(s,h,p) 


(52) 


which  after  inversion  yields  the  behavior 


u(x,  h,  p) 


(53) 


)  =  0(x2) 

v(x,  h,  p) 

The  actual  terms  corresponding  to  (52)  [and  (53)]  are  easily  derived. 
Without  the  singular  terms  in  (47)  the  remaining  terms  there  behave  as 
(l-y/h)x  near  y=h.  So  if  (1-y/h)  gets  small  at  the  same  rate  as  x  does,  this 
behavior  is  in  agreement  with  that  in  (53).  The  singular  terms  in  (47), 
of  course,  get  large  as  y-*h  which  is  not  consistent  with  their  absence  at 
y=h.  This  suggests  that  additional  terms  of  the  singular  type  in  (47)  might 
be  involved,  all  of  which  sum  to  give  zero.  Indeed  carrying  the  parts 
integration  of  u^  on  further,  for  example  (: v ^  is  similar),  shows  that  for 
small  x  it  has  the  singular  terms 


upi(x,y,p) 


a0<P) 
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(54) 


After  the  first  term  in  (54)  there  is  sign  alternation  of  the  terms  in  each  of 
the  brackets,  and  this  pattern  is  true  for  the  higher  order  terms  also. 

TT 


The  order  term  in  (45)  is  not  present  in  the  y=h  solution. 
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Hence,  it  is  reasonable  to  expect,  that  as  y-*h,  the  sum  of  the  terms  in 

(54)  will  approach  zero.  It  also  follows  then,  for  the  leading  singular 

term  in  (54)  to  represent  u^  in  (47),  as  it  does,  x  must  vanish  at  a  faster 
1/4 

rate  than  (1-y/h)  does,  i.  e.  the  near  field  solution  (47)  is  limited  to 
smaller  and  smaller  x  values  as  the  corner  at  y=h  is  approached.  Other¬ 
wise  the  use  of  the  additional  terms  in  (54)  will  be  needed,  which  obviously 
will  still  have  limitations  as  y  gets  closer  and  closer  to  h.  One  can,  of 
course,  evaluate  at  a  station  quite  near  y=h  in  this  manner,  then  evaluate 
at  y=h  from  the  terms  of  which  give  (53). 

It  is  clear  then,  that  once  the  unknowns  a^p),  an(p)>  and  b  (p) 
are  numerically  evaluated,  a  simple  inversion  of  (47)  [since  it  has  already 
been  shown  the  inverses  of  the  unknowns  behave  as  H(t)]  will  yield  the 
long  time-near  field  solution  for  the  present  problem. 


Verification  of  Long-Time  Solution. 

Short  of  writing  the  inverse  (26),  which  will  be  postponed  until  the 
numerical  work  for  the  coefficients  of  the  edge  unknowns  is  carried  out 
(as  discussed  earlier),  some  things  can  be  done  on  verifying  that  (15), 
defined  by  (35)  and  (36),  will  yield  a  solution  satisfying  the  present  boundary 
value  problem,  i.  e.  the  differential  equations  (1),  the  initial  conditions  (3), 
the  boundary  conditions  at  the  edge  x=0,  (21),  the  plate  face  y=h,  and 
center  y-0,  the  inverses  of  (10)  and  (11),  and  the  conditions  at  x-*oo,  (4). 

Verification  that  the  long-time  solution,  once  it  is  written,  satisfies 
the  differential  equations  (1)  can  be  shown  by  direct  substitution.  Since 
the  problem  is  based  on  (1),  a  system  of  hyperbolic  equations,  one  would 
expect  the  general  solution  (15),  defined  by  (12)  and  (8),  to  have  a  singularity 
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free  half  plane  Repa  Y,  a  positive  constant.  Hence  completion  of  Br^  in 
this  domain  would  give  the  corresponding  solution  a6 


u(x,  y,  t) 


v(x, y,  t) 


=  0  ,  for  0<t<  Vx  +y  / c 


hence  satisfying  not  only  the  initial  conditions  (3)  on  the  displacements, 
but  also  on  the  velocities  there. 

Direct  verification  of  the  edge  boundary  conditions  (21)  can  be 
obtained  by  setting  x=0  in  (47)  and  (53).  A  differentiation  of  these  w.  r.  t.  x, 
and  then  letting  x-*0  [noting  the  order  terms  in  (44)  to  (46),  and  the  higher 
order  x  terms  in  (54)1,  recovers  the  assumed  expressions  for  the  strains 
u  (0,  y,  p)  and  v  (0,  y,  p)  in  (22).  Note,  in  the  case  of  (54),  after  this 

X  X 

differentiation,  that  setting  x=0  away  from  y=h,  leaves  only  the  leading 
term  as  it  should  be.  But  when  y-*h,  it  is  now  reasonable  to  expect  that 
the  terms  only  after  the  first  (which  in  each  of  the  brackets  alternate) 
would  sum  to  zero. 

Direct  substitution  of  y  =  0  in  the  second  of  (13),  with  the  second  of 

(35)  and  (36)  defining  the  integrand  there,  gives  v(x,0,p)=0,  hence  from 

(15),  v(x,0,t)  =  0  (independent  of  time,  hence  giving  the  same  result  from 

(26)  for  long  time).  Thus  the  double  inverse  of  the  second  of  (11)  is 

satisfied.  Now,  using  (13),  with  (35)  and  (36),  to  form  o  (x,y,p)/ti,  in 

xy 

accord  with  the  third  of  (2),  it  is  easily  shown  by  direct  substitution  of 

y=0  that  a  (x,  0,  p)  and  therefore  a  (x,  0,  t)  are  zero  Thus  the  double 
7  xy  xy 

inverse  of  the  first  of  (11)  is  also  satisfied. 

Now  substituting  y=h  directly  into  the  Br  integral  just  formed  for 

s 

o^fx,  y,  p)/(j  ,  and  a  similarly  formed  one  for  a^x,  y,  p)/ (X+2^  ,  in  accord 
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with  the  second  of  (2),  it  is  easily  shown  that  a  (x,h,  p)/|i=0  ,  and 

xy 

Oy(x,h,p)/(X+2u)=a06(x-a)G{p)/(X+2^)  ,  and  hence  a^(x,  h,  t)/n=0  ,  and 
Oy(x,  h,  t)=cr06(x-a)G(t).  Thus  the  double  inverses  of  (10)  [see  also  (20)a] 
are  satisfied. 

It  remains  to  show  the  conditions  at  x-*oo  ,  (4),  are  satisfied.  In 
the  case  of  the  general  solution  (15),  with  (35)  and  (36),  again  (55)  assures 
satisfaction  of  (4)  in  advance  of  the  wave  front,  i.  e.  for  /x  2+y2  >cdt. 
However,  back  of  this  wave  front,  and  for  the  harmonic  problem  solution, 
(13),  with  (35)  and  (36),  it  is  necessary  to  show  these  solutions  do  not 
become  unbounded.  The  Tauberian  theorem,  in  which  s-0,  can  be  used 
to  show  that  this  far  field  condition  is  satisfied.  There  are  two  cases  to 
be  considered.  First,  s-*0  and  p  small  but  not  zero,  and  second,  s  and 
p-*0  together.  Approximations  of  (35)  and  (36)  corresponding  to  the  first 
case  show,  at  most,  contributions  of  0(1)  in  s  for  0  ^  y  £  h.  It  follows  that 
u(x,  y,  t)  and  v(x,  y,  t)  behave,  at  most,  as  6(x)  for  x  -»  oo,  and  therefore 

C3M  2* 

vanish.  Since  u^,  v^,  and  higher  order  derivatives  in  x,  involve  higher 
order  behaviors  in  s  (as  s-»0),  it  follows  their  inverses  are  also  zero  for 
x-*ao.  u  ,  v  ,  and  higher  order  derivatives  in  y,  differ  from  u  and  v 
only  through  higher  powers  of  a  and  (3,  and  since  for  s-*0  the  latter  behave 
as  k^  and  kg,  respectively,  the  inverses  of  these  derivatives  behave  like 
u  and  v,  vanishing  as  x-*oo.  It  follows  that  for  the  present  case  (4)  is 
satisfied  everywhere  for  large  x. 

Now  for  the  case  of  s  and  p-*0  together,  R  behaves  as 


R  =-  -b2hpkg2  [s2-kp2] 


(56) 
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222222  /  j  ' 

where  b  =c^  /cg  '  =P  ^Cp  *  Cp  being  the  "plate  velocity"  /E/p(l-v‘:)  . 
(56)  is  well  known,  the  zeros  of  ^s^-kp^J  therein  corresponding  to  the  low 
frequency-long  wave  approximation  to  the  lowest  mode  of  R  in  (14),  i.  e.  , 
the  plane  stress  approximation.  Approximating  (35)  and  (36)  for  this  case, 
shows  that  for  u(s,v,p) 


R.  . 

uc  i  (s  ’  y,  p) 


b“  s  -k 


n-1 

)  (-f  -  a  (p) 

Zj  mr  n  r 


p  /  n=l,3,5,.. 


uc2(s,y,p) 


k  XaQ(p) 
hbVs2-k 


survive  for  O^y^h,  where  X-C^  =Ck  .  Now  since  (33)  shows  the  coef- 

d  Ks 

ficients  aQ(p)  and  afi(p)  [and  bn(p)l  behave  as  1/p,  for  G(t)=H(t)  the  step 
input,  simple  table  inversions  of  (13)  and  (15),  with  ua“uci  +  uc2  where 
the  latter  are  given  by  (57)  and  (58),  yield  the  behavior  of  u(x,y,t)  to  be 

u(x,  y,  t)  =-  U(t  -  x/Cp  )  (5 

15 

for  the  long  time-far  field  solution  ,  where  U  is  a  constant. 

Similarly,  for  v(s,y,p)  in  0<y^h,  approximation  shows  that 


vcl(s.y.p): 


(2-k  j  ys 


n-1 

T  2 


b  ( s  -  k 


Y  (-f  —  a  (p) 

L,  nir  nvr/ 


p  /  n=l,3,5,.. 


vc2(s,y,p)  = 


(2-k^)y  sXaQ(p) 


b"  s“-k 


survive,  which  leads  to  the  double  inverse 


v(x ,  y,t)  -  Vy  H(t-x/cD  ) 


An  equal  term,  but  corresponding  to  a  negative  traveling  wave,  is  ruled 
out  in  the  inversion  of  (15),  since  it  leads  to  an  unbounded  contribution. 


for  the  long  time-far  field  solution,  where  V  is  a  constant.  It  should  be 
noted  that  (59)  and  (62)  are  restricted  to  a  small  region  just  after  t=x/cp. 
There,  however,  they  correspond  to  non-decaying,  in  space  (x)  and  time, 
bounded  disturbances.  Because  of  the  latter  these  solutions  satisfy  the 
requirement  of  boundedness  in/x  +y  cc^t  ,  since  CpCc^.  The  non¬ 
decaying  nature  of  this  disturbance  has  been  observed  before,  in  several 
related  problems,  e.  g.  the  work  of  Curtis  [  1 5 1 .  It  is  of  interest  to  note 
that  the  time  behaviors  ofuandv  in  (59)  and  (62),  at  a  plate  station,  agree 
with  the  time  behaviors  of  the  corresponding  radial  and  thickness  dis¬ 
placements  in  the  solution,  given  by  the  author  [16],  for  the  related 
problem  of  excitation  of  an  elastic  plate  by  symmetric  normal  point  loads 
(see  Figs.  4,  5  in  [  16]  ).  The  latter  displacements,  however,  decay 
spatially.  It  is  planned  to  include  a  more  complete  Btudy  of  the  lowest 
mode  disturbance  in  the  future  work  on  the  present  problem. 

§4.  COMMENTS  AND  FUTURE  WORK. 

Plans  to  evaluate  further  the  solution  treated  here,  have  already  been 
discussed.  It  is  of,  at  least,  equal  importance  to  remind  the  reader  that 
solution  of  the  integral  equations  (19),  basic  to  this  work,  has  so  far  been 
limited  to  the  eigenfunction  technique  employed  by  Benthem  in  his  problems 
from  elastostatics  [see  (22)].  Their  applicability  here  has  been  argued  on 
the  basis  of  the  long  time-low  frequency  analog.  It  is  therefore  desirable  to 
look  into  other  possible  more  general  means  of  evaluating  (19),  so  that,  per¬ 
haps  short  time-high  frequency  phenomena  can  be  investigated  in  the  problem 


treated  here. 
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Further,  solution  of  (19)  for  problems  of  the  other  non-mixed 
edge  type,  i.  e.  free  edge,  including  sudden  edge  loads,  is  also  of  strong 
interest.  Some  work  has  been  done  on  these  problems  and  plans  are  to 
continue  it. 
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